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Abstract 

A classic sum rule by Das et al. is extended to seven of the low-energy constants Ki, 
introduced by Urech, which parameterize electromagnetic corrections at chiral order 0{e^p^). 
Using the spurion formalism, a simple convolution representation is shown to hold and the 
structure in terms of the chiral renormalization scale, the QCD renormalization scale and 
the QED gauge parameter is displayed. The role of the resonances is studied as providing 
rational interpolants to relevant QCD n-point functions in the euclidian domain. A variety of 
asymptotic constraints must be implemented which have phenomenological consequences. A 
current assumption concerning the dominance of the lowest-lying resonances is shown clearly 
to fail in some cases. 
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1. Introduction: 


Thirty years ago, Das et al.[|I| (DGMLY) derived a remarkable relation between the mass 
difference of the charged and neutral pions and the masses of the lightest vector and axial vector 


resonances, 


m2 - = 


3e2 


MlM^r 


167r2 Ml - Ml 


log 


Ml 

Ml 


( 1 ) 


This relation follows from a sum rule which is exact in the chiral limit (i.e. = rud = = 0) 

under the only extra assumption that the lowest-lying vector and axial vector meson resonances 
make the essential contribution to the integral. The physical mass difference is nearly 

purely electromagnetic in origin and happens to be rather accurately described by eq.([^). The 
analogous mass difference of kaons, M^+ — Mj^o, has an electromagnetic contribution and 
a purely QCD contribution proportional to — rrid which are approximately of the same 
magnitude. Knowledge of the electromagnetic contribution allows one to access the value of 
the quark mass difference rriu—md (divided by, say, using chiral perturbation theory^. 

It has long been believed that estimating the EM contribution to M|.+ — M^-o in the chiral limit 
was sufficient. In this limit, it is given by Dashen’s theorem P] (DT) to be equal to Ml+ — M^q. 
This approximation is now known to fail for the purpose of extracting — rrid. In particular, 
the value of the rj decay rate Tij] ^ Svr) that one would predict (at one loop and including 
estimates of higher loop corrections [Q|Q) would be too small by as much as a factor of two 
compared to experiment. 

Thus, it is necessary to estimate the EM contributions to M|.+ — Mj^o beyond the chiral 
limit (the difference from Ml+ — M^o is customarily, but somewhat inappropriately, referred to 
as the violation of Dashen’s theorem). There have been many attempts over the years in this 
directionp0 (a representative, but non-exhaustive list) and, very recently, a 

lattice calculation has appeared|^. Several rather different approaches to this problem have 
been followed. The work of refs.0 and [O is based on the assumption that most of the 


information on DT violation is contained in the parameters of the low lying resonances, in a 
way similar to eq.(|lD. Unfortunately, the results of these two papers (which should be identical) 
are not in very good agreement which each otherQ. One of the purposes of the present paper is 
to investigate in some detail the validity of this assumption. 


As shown by UrechfT^, at low energies, electromagnetic effects can be parametrized 


m a 


chiral lagrangian framework |T^ [§]. In this formalism, the leading electromagnetic contribu¬ 

tions are computed from the effective lagrangian at tree level, which contains a single low-energy 
constant (LEG): the result of Das et ah provides a sum rule for this LEG. At next-to-leading 
order, the electromagnetic corrections are obtained by computing the photon loop (as well as 
the pion loops) from the leading order lagrangian and adding the contributions at tree level 
from the next order terms in the lagrangian, which involves (essentially^ ) 13 new low-energy 
constants Ki, K 2 ,Ki^. In this paper, we will propose a sum rule evaluation of the seven 
parameters Kj ,..., K 13 . The mass difference actually involves two more constants, 

and iLg, which we will not attempt to evaluate. Little was known, up to now, on the individual 


2 A diagram was incorrectly evaluated in ref.Q. An update of this calculation has appeared very recently pS 
which obtains a very similar result numerically. 

^There are further constants which either play no role at low energy or correspond to contributions. 
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values of these LEC’s (estimates for some combinations were given recently [T^) except for their 
order of magnitude which, for consistency of the chiral counting for the electric charge, should 
be the same as that of the ordinary O(p^) constants Lj. Individual knowledge of these constants 
is useful for the purpose of evaluating radiative corrections. Indeed, as precision increases in the 
predictions of chiral perturbation theory, the computation of radiative corrections is becoming 
a topic of increased interest. This is already true for the tttt scattering amplitude which is now 
known at the level of two loops|^|^, and which will be subjected to a precision test in the 


planned experiment DIRAC [IS 


In general, we will show that the iCj’s can be expressed as a convolution of a QCD correlation 
function with the electromagnetic propagator, plus a contribution from the QED counterterms 
which remove the divergence of the integral. The LEC’s Ki,...,Kq are related to QCD 4- 
point functions, while Kj ,..., iLia are related to QCD two- and three-point functions. In that 
case, we will show that the contribution of the resonances can be discussed independently of 
any specihc lagrangian model for resonances. An important question concerns the validity of 
the approximation of retaining only the lightest multiplet of resonances in each channel. One 
may view resonance saturation as a method of constructing rational interpolants to the QCD 
n-point functions in the euclidian region. These interpolants must satisfy certain asymptotic 
constraints, in order for the QED divergencies to cancel out, and one also expects them to be 
reasonably precise at low momenta, i.e. to some extent in the resonance region as well. It is 
not clear that low order interpolants are capable of satisfying all these constraints. 


The plan of the paper is as follows. In the next section, we discuss the form of the QED 
counterterms when the spurion formalism is applied to the electric charge. The use of this 
formalism is a key technical ingredient in the derivation of the sum rules. In sec.3, we discuss a 
set of sum rules for the parameters Kn, K 12 and iCia which involve QCD correlators in the chiral 
limit. We will show that a single multiplet of vector and of axial-vector resonances is enough 
to obey adequately all the asymptotic constraints, and we obtain expressions which are rather 
neat generalizations of eq.(|^). Next, in sec.4, we discuss Kj, Kg, Kg and Kiq. The hrst two are 
suppressed by the Zweig rule and the latter two can be expressed in terms of flavour symmetry 
breaking differences of vector and axial-vector spectral functions. The phenomenology in terms 
of resonances is then discussed in sec.5 and the application to the violation of Dashen’s theorem 
in sec.6. 


2. The spurion formalism and the QED counterterms 

The spurion formalism applied to the electric charge matrix |]^ q = diag(2/3, —1/3, —1/3) 
consists in calling q by two different names qi and q^ such that the QED interaction lagrangian 
for the quarks can be split into two pieces: 


+ iq4)'^ = + iqL4)i^L + + iqR4)^R ( 2 ) 

(with 4’^ = {u,d,s) ). One may thus render the QED lagrangian invariant under the chiral 
group by assuming appropriate chiral transformation properties for the two spurions qi and 
qR. This is the same method as has been applied to the mass matrix|]T^. In order to take 
quantum electrodynamical effects into account, the chiral lagrangian must include the photon 
as a dynamical held, in addition to the octet of pseudo-Goldstone boson helds. As usual, 
the chiral lagrangian will consist of the most general set of local interaction terms which are 
invariant under the chiral group and which are classihed according to increasing chiral order. 
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In the EM sector, the correct counting rules, introduced by Urech[14], state that the photon 


held Afj, is of order 0(jP) and the charge spurious qi and qu are of order 0{p). At order two, 
the chiral lagrangian consists, hrstly, of the ususal terms. 


£(2) = ^ < D.UD^U^ + t/yt + 


(3) 


where, however, the covariant derivatives contain qi and qR in addition to the usual vector and 
axial-vector external sources 


Df^U = d^U - i{v^ + a^ + qRA^)U + iU{v^ -a^ + qiA^). (4) 


In addition to that, one must include the purely photonic lagrangian. 


p(2) 


_ J? 

’ 


and, hnally, there is a single invariant term with two spurious that one can form[pO 


=C< qRUqLU^ > 


(5) 

( 6 ) 


Equations (i, (D, (I) dehne the most general lagrangian of chiral order two which satsishes 
EM gauge invariance in the limit of constant spurious. In order to quantize the electromagnetic 
held, one must include a gauge hxing term and it is also useful to introduce a small photon 
mass My in order to regulate infrared divergencies at intermediate steps of the calculations, 
this is performed by adding 


’(2) _ ^ ]\/r2 - - 1 


=-M^ A,A^ - - {d,A^y . 


The corresponding photon propagator is 

-iD^yp) = 


pZ 


M2 


2^ 


dfiu + ('C “ 1) 


PflPv 


- eM2 


(7) 

( 8 ) 


At chiral order two. Green’s functions and observables are computed by using at tree level. 
At the next chiral order, O(p^), one must hrstly compute the one-loop graphs with vertices 
generated from The loops will involve both pion and photon internal lines in the present 
case. Secondly, one must add the tree contributions from The lagrangian will consist, 
hrstly, of the usual terms of Gasser and Leutwyler.]^]. In addition, one can construct a set of 
terms with two EM spurious 0, which were classihed by Urech. For convenience we will separate 
the terms which involve two derivatives of the pion helds, 

X q^qL + > (9) 

+A2F2 < D^UD^U^ X qLU^qRU > 

+K^F^{< D^UqRU^ >2 + < D^UUrU >y 

+A4F2 < D^UqLU^ X D^U^qRU > 

+K,F^ < D^U^D^UqRqL + D^UD^U^RqR > 

_ +K^F^ < D^U ^D^UqiU^qRU + DqRUqRU^ >, 

^Symmetry considerations alone do not rule out terms with a single spurion and no for example a term 
like < qLDfj,WD^U + qRO^UD^^V^ >. The technique to be used repreatedly below immediately shows that 
the coefficients of such terms vanish because the vacuum expectation value of the photon field does. 
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from the remaining ones 


< X^U + U^x >< QlQl + qRQR > ( 10 ) 

+K^F^ < x^U + U^x >< qiU^qRU > 

+KgF^ < {x^U + U^x)qLqL + {xU^ + Ux^)qRqR > 

+K,oF^ < ix^U + U^x)qLU^qRU + {xU^ + Ux^)qRUqLU^ > 

+KuF^ < ix^U - U^x)qLUkRU + ixU^ - Ux^)qRUqLU^ > 

+KuF^ < UDf^U^[D^qR,qR] + D^U[D.qR^qj] > 

+K^^F^ KD^qRUD^^qLU^ >, 


where 

D^^qi = d^qi - i[Vf, - a^, ^l], = d^qR - i[Vf, + a^, qR] . (11) 

Fnrther terms which contain fonr EM spnrions as well as terms which vanish when is set to 
zero will not be considered in the following. 


All these terms are connted as O(p^) according to the standard chiral connting rules ||I6[|. 
According to the so-called generalized CHPT (see e.g. ref.[^) quark masses are counted as 
0 {p) rather than 0 {p^) such that the terms corresponding to iPy up to Ku would be counted 
as of chiral order three. In this scheme, the terms corresponding to Ki to Kq as well as K 12 , 
Ki 3 would appear at O(p^) as before but, to the same order, a number of additional terms 
with two occurences of the scalar held x must also be included. This modihcation of the chiral 
counting corresponds to a weak form of spontaneous chiral symmetry breaking in which F^^ ^ 0 
while the quark condensate < qq > could be small or even vanishing. Whether this scheme 
is of physical relevance or not can be decided experimentally by studying the low-energy tttt 
amplitude, which is very sensitive to the size of the quark condensate p^|^. In this paper. 


the discussion will be restricted to the case of the standard chiral counting, which is the most 
predictive framework. 


The effective lagrangian generates a low-energy representation for the generating functional 
of QCD Green’s functions W. In addition to the sources s{x), p{x), v^{x), a^{x) that one 
customarily introduces, W now depends on the two new sources and qR, 


fh = W{v^,a^,s,p,qL,qR) ■ 


Accordingly, one can introduce Q-currents, 


Q 


a _ 

L — 


6 W 

iSql' 


Qr = 


6 W 
i5q% ■ 


( 12 ) 


(13) 


and one may identify the low-energy constants by taking functional derivatives of W with 
respect to g^ and qR. It is useful, in connection with the EEC’s 7^12 and iPis, to further extend 
this formalism by letting g^ and qR become space-time dependent. This generalization allows 
one to identify iEis, for instance, from the two-point function 5‘^W/5qL{x)5qR{Q) instead of a 
four-point function if the spurious are restricted to be constants. There is an apparent drawback 
to this extension, which is that we are allowing the photon to couple to a nonconserved current. 
This is an unusual situation and such a theory is not guaranteed to be renormalizable. It is not 
clear whether it is possible to dehne the generating functional in full generality in the presence 
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of X dependent spurions. For our restricted purposes (i.e. the derivation of sum rules which 
are UV hnite for K 12 and iFia) we simply need the expansion of W to quadratic order in qi, 
qn and we can switch off the external electromagnetic held. In this case, the UV divergencies 
are removed by the following hnite set of counterterms 


^^ _ _ 

-^92 = -Z2i)LqLJ^qL'4^L + {L^ R)- Zsi}RqR{s + ip)qL'4^L - Z.'ilJiqLis - ip)qRipR (14) 


The existence of such a set of local counterterms guarantees the validity of the idea of using 
space-time dependent spurions at this order. For our purposes, (|^) will be used at tree level. 
We may then use the equations of motion generated by (|^), thereby obtaining the counterterms 
in the form 


D^qL]-ff,'ipL + {L^ R) 


(15) 




-^2 

2 


s zp) -f (s ip)ql) - ZsqR{s ip)qL 


^^2 (qR{s - ip) + (s - ip)ql) - Z^qLis - ip)qR 


'^L 

'tpR ■ 


Setting qi = qR = q, one recovers the usual QED counterterms. For a free fermion, the 
renormalization constants Z 2 and Zg are determined from the requirement that the fermion self 
energy F(ji) is normalized in the way appropriate for a free particle of mass m. For instance, 
in dimensional regularization, Z 2 and Zg read 


, -2e f 2 'j 

I -dr(') +1 - iy;|) + K - 3) + Oog{ - sj (w) 

, -2e ( 2 1\X2 'j 

+ +1 - log + (e - 3) log ^ + eloge + 2|, 


Quarks, however, are not free particles but on the contrary, are conhned, so it seems unwise to 
use the above values of Z 2 and Zg in the present context. In fact, the QED counterterms play 
a role in a kinematical region where QCD becomes perturbative. The most reasonable thing to 
do, then, seems to resort to MS renormalization, as one customarily does in perturbative QCD 
calculations. Such a mass independent renormalization scheme is in fact implicitly assumed in 
the dehnition of the “current” quark masses, which appear in the effective lagrangian (see e.g. 
the review by Manohar in the PDG[^). As a consequence, the corresponding renormalization 
scale /io will appear in the expression of some of the low-energy constants. This was pointed 
out in ref. |jl^ . According to this prescription, Z 2 and Zg have the following expression. 


Z'2 = {^(-?)[r(£) + log(47r) 
" levr^' 


- 3) F(e) -F log( 47 r) 


(17) 


which will be adopted in the following. Further divergences involving the strong coupling 
constant ag will be assumed to be removed by the MS prescription as well. 
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3. QCD in the chiral limit: sum rules for ( 7 , K 12 and Ki^, 

3.1 The constants C and spectral representation 

Let us consider the correlation function < QlQr > constructed from the currents associated 
with the spurion helds qi and qR introduced in sec.2, 

n?H(p") j <(,\TQl{x)Q\m(>> (18) 

= i Jd‘xe‘f‘ < 0\TQl.{x)Ql,{0) - Q5,(i)Q;i(0)|0 > 

We have introduced vector and axial vector Q-currents: 

Qv = Ql + Qr Qa = —Ql + Qr ■ (19) 


On the one hand, we can use the chiral lagrangian, and compute the chiral expansion of 
Ulr up lu chiral order four. It receives contributions at tree level involving the constants C and 
iLi 3 and a photon loop contribution. There is no pion loop contribution. Explicitly, a simple 
calculation gives, in an arbitrary gauge^ 


n«s(p 2 ) = 2 C + pifi 


2KUp) + 


1 3 


Ibvr^ 4 


(?-!)( 


, 1 

log ^ + - 


+ '^log^ 


+ 0{p^ 


( 20 ) 


On the other hand, it is clear from (H) that each Q-current is nothing but a QCD current 
multiplied by the quantum photon held, so that the above correlator can be written as a 
convolution of a QCD correlation function with the photon propagator: 

n?R(p") = y {K'iik - p) - n“?(fc - p)) . (21) 

In this case, the QED counterterms (see (|l^)) make no contribution, implying that the integral 
must be hnite. This can be verihed by using the operator product expansion for the difference 
n ^3 — which, indeed, starts with dimension six operators in the chiral limit. The vector 
and axial-vector currents are dehned, as usual, as 


yixix) = Al{x) = . ( 22 ) 

The < VV > and < AA > correlation functions which appear in (^) satisfy once-subtracted 
dispersive representations (which we display in unsubtracted form for simplicity and for the 
reason that we will only consider combinations which need no subtractions in the sequel). 
Assuming the T-product to be covariantly dehned, 

/I fOO r] o 

d^xe*^^ < 0|Tl/3'^(a:)l/370)|0 >= (gV -d'^V)- / - o Pv^s) (23) 

TV Jo s — q^ 

and 


I{%{q) = i J d^xe*"^ < 0\T A^>^{x)A^’^{0)\0 > 

= {q^q"'- g^''q‘^)- f ^ pA^i^) + 9^’'- j ^^c^A3(s). 

IT Jo s — q^ IT Jo s — q^ 


(24) 


^ We follow the usual chiral MS convention of relating bare and renormalized constants by: Ki = 
- Il7327r2(r(e) + 1 + log(47r))] 
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In the chiral limit, the spectral function CAsis) vanishes identically. It is convenient to separate 
explicitly the pion contribution in the axial spectral function, 


Pa»(s) = Pa>(s) + . 


(25) 


It is now a simple matter to insert the spectral representations (|^) into the expression (|^) 
for and compute the photon loop integral. Comparing with the chiral expansion form of 
nfij (0)) 0116 readily identihes the constants C and Ki^ for which one obtains the following 
represent at ions: 

2 roo ^ 

C = dxx\og—{pv3{x)-pA3{x)), 

'Jo 


(26) 


27rJo ^ p^ 

which is a re-derivation of the DGMLY sum rule[0] (using the relation M^+ — M^o = 2e^C/F q 
which is exact in the chiral limit and to order 0(e^)), and its immediate generalization. 





?) 


1 

12 ^ 


pOO ^ 

/ dx log ^ {pv3 (x) - pa3 (x) ) 
Jo 


(27) 


The expected order of magnitude for the constants Kl{p = My) ~ is indeed conhrmed 

by this explicit calculation, but one must be careful since these constants are gauge dependent 
in general. The dependence of C and Ki^ upon the chiral renormalization scale p can be 
deduced from eq. (^OD since the left-hand side is independent of p\ C is a constant and iCis 
satishes 




dKM 


1 3 


K -1). 


(28) 

and ( 0 ) provided 

the spectral functions satisfy the two Weinberg sum rulesp^ (a modern discussion can be found 
in ref, ) 


dp IGvr^ 4 

The correct scale dependence emerges from the spectral representations 


POO POO 

— dx {pv3{x) — pa3{x)) = Fq, / dxx{pv3{x) — pa3{x)) = 0 . 

TT JO Jo 


(29) 


3.2 The EEC’s Kn and K 12 and the correlator < VAP > 

We will follow the same approach as in sec. 3.1. In order to isolate contributions from 
the EEC’s Kii and K 12 we may consider a correlation function of two Q-currents with the 
pseudoscalar current, 

n^(p,g) = / d^xdSpp^+^^y < 0\TQUx)Q\{y)P%0)\0 > , (30) 

where the pseudoscalar current is dehned as 

P°‘{x) = iV'(x) 7 ^^-^(x) . (31) 

It will also be useful to introduce the related pion to vacuum matrix element, 

ng.feO = / d‘xe‘'^ < a\TQ\,{x)Q\(y)\AHl) > 
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(32) 

















We first construct the chiral expansions of these objects, up to one loop, using the chiral 
lagrangian. Setting g = 0 in (|50[), we obtain 


n^(p, 0 ) — 


-2BoC 

p2 


■ 2BoF^ {2Kn - + 0{p^) . 


(33) 


The contributions from the photon loop and from the pion loop both vanish, implying that 
the combination 2Kii — is hnite in any gauge. Setting p = 0 now in (^0|), we obtain the 
following chiral expansion, which involves the combination 2Kii + i^i 2 . 


nS(0.9) 

qz 


(34) 


-BoF^ h{2K[,(ii.) + K[, 

We will also use the expansion of 11^: 


+ 




3 1 


+ 0{ct 


2C ( 

nO (p,0 iFo.P ■' { - 4A'[j(p) - 2KIM 


0 

1 1 


(36) 


Ibvr^ 4 


1 Q-i > 

(-7^ + 9) log ^ “ 2^ ^ 2 ^ P ■ 0| + 0{p^) 


We may now express 11^ and 11^ by functionally differentiating the QCD generating functional 
with respect to the sources Ql^x), Qr^x) and p{x), and one obtains. 


r d'^k 

n%{p,q) = J j:^,W^.{p-k,q + k)(-,D'‘-'(k)) + 




Z2~\- Zg — 2Z' 


p-{p + <l) 

(p + qV 


, (36) 


where W^'^{p,q) is the correlator of one vector, one axial-vector and one pseudoscalar current 
< VAP > in the chiral limit, 

W^,{p,q) = J < 0\TV^{x)Al{y)P\0)\0 > . (37) 


In eq.(^) the convolution integral simply arises from replacing the currents Qv and Qa from 
their dehnitions. The additional contribution gets generated upon functionally differentiating 
the QED counterterms. This contribution implies that the convolution integral is UV divergent 
in QCD, in such a way that the sum of the two terms in eq.(p6D is hnite. 

Equating expression (|3^) with the chiral expansions (^3]) and (^4]) one obtains an expression 
for the EEC’s Ku and Ku as a convolution of the QCD 3-point function < VAP > and the 
electromagnetic propagator. An additional expression, which we have used to cross-check the 
results, can be found which involves the pion to vacuum matrix element < 0|I7A|7r >, 


W^ApJ) = f d^xpP^ < 0\TV^{x)Al{0)\7r^{l) >= —^\im fW^^{p,q) . (38) 

The Q-current matrix element introduced in eq.(^) can be expressed in terms of < 0|I7A|7r > 
as 

c /r ^ 

(p, l) = J 7^M>(p - k, in-W^^k)) + iFqZ^p-1 , (39) 
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and using the chiral expansion (^) one obtains a convolution representation for the combination 

2Ki2 + i^l3. 

To summarize, at this point we have shown that the constants i^n, K 12 and Ki^ can 
be expressed as convolution representations in terms of QCD correlators in the chiral limit. 
Using the same ideas it is not difficult to see that the constants Ki to Kq must obey similar 
representations, involving the QCD correlators < VVAA > and < AAAA > in the chiral limit 
as well. The difficulty in this case is the implementation of resonance saturation. We will now 
show that a rather neat implementation can be performed for the correlators VV — AA (which 
is very simple) and < VAP > ( which is somewhat less trivial) and defer a brief discussion of 
Ki to Kq to sec. 3.6. 

3.3 Resonance saturation of < VV — AA > 


The spectral densities pv and pA which appear in the expressions (|2B|) and (P?]) for C and 
/C 13 can be related to exprimentally measurable quantities like > hadrons cross sections 

and r decay rates. Analyses of the DGMLY and Weinberg sum rules in terms of available 
experimental data have been performed [^[^. A physically appealling and time honoured 
approximation is that of “resonance saturation” of the spectral integrals. There are, in fact, two 
aspects in this approximation. One is to represent the spectral functions as a sum over delta 
functions, and the second is to retain just a few terms in the sum. The hrst aspect corresponds 
to a leading large limit and was actually shown to be rather precise in practice However, 
even in the large Nc limit it is not clear how many resonances should be included. A current 
assumption is that minimal saturation, i.e. in the present case including only the contribution 
of the p(770) and the ai(1260) resonances, does provide a resonable approximation. In this 
case, the spectral functions are given by 


PV,M = Fi-S(x - Mp, Px>(x) = FlS(x - Ml) 


(40) 


One may at least verify that this lowest order approximation is not grossly inconsistent. Indeed, 
if one uses experimental values for Fv and My then the Weinberg sum rules are satished 
with values of Fa and Ma which are within 30% of those of experiment. Furthermore, the 
corresponding result for the constant C seems to be rather accurate as one can judge from the 
resulting value of M^+ — M^o • In the dispersive representation, one may attribute the success 
of the minimal resonance approximation to the fact that the various integrands are oscillatory. 
In that situation, numerical analysis teaches us that an optimal sequence of approximations is 
obtained by cutting the integrand after an even number of oscillations. The precision could 
be much better than one could expect from the asymptotic behaviour of the absolute value of 
the integrand. This also suggests that in order to improve on the hrst order approximation 
one should include an additional pair of resonances p' and aj, which means integrating up to a 
rather high energy. 

An alternative point of view on resonance saturation emerges from considering the corre¬ 
lation function < VV — AA > itself. Resonance saturation amounts to constructing rational 
parametrizations for the correlation function. In the minimal version, for instance, one would 
have 


n^ys(p)-n^A‘'4p) = (p^p‘'-9'^y^ 


p2 


+ 


Fl 


P2 


p'^-Ml 


p'^-Ml 


+ polynomial 


(41) 
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This is the most general form which satisfies the chiral Ward identities (which dictate the 
form of the tensor strnctnre) and which has a specified nnmber of resonance poles. When 
computing the four-dimensional convolution integrals for C or it'is, one may perform a Wick 
rotation such that the correlator < VV — AA > is needed only for euclidian values of p^. In 
this region, the correlator is a smooth function so it makes perfectly good sense to employ a 
rational approximation. The parameters of the rational function (|4l|) (My, Ma, Fy,...) may be 


constrained both from the asymptotic region > — cx), using information from the operator 
product expansion, and by extrapolating to positive values of p^ using physical information 
on the resonances. In the asymptotic region, the behaviour must ensure that the convolution 
integral defining 11®^ eq. (|^) is finite. This implies that the polynomial in eq.(|4l|) must vanish 
and one must impose the two conditions that the coefficients of 1/p^ and 1/p^ both vanish. 
These conditions are exactly equivalent to the two Weinberg sum rules. We recover here 
the necessity of imposing asymptotic matching conditions when discussing resonances in the 
effective lagrangian context, a fact which has been appreciated relatively recently 0 (see also 
pop. The usefulness of the effective lagrangian for the resonances is to guarantee that the 
n-point functions that one might compute automatically obey all the chiral Ward identities. 
When discussing two- or three-point functions it is sometimes more expedient to write directly 
the most general form compatible with a given set of resonance poles and impose that the 
relevant Ward identities be obeyed by hand. This is what we have done for < VV — AA > and 
we will proceed in this way also for < VAP > below. 

Returning to the evaluation of the low-energy constants, after imposing the two asymptotic 
conditions on the parametrization (^T]) , performing the Wick rotation and integrating over 
angles one obtains C in the form of an integral 


^7 = 




dx 


327r2 


{x + M^){x + Ml) 


(42) 


The interest of this expression is the observation that a substantial part of the integral (ap¬ 
proximately 40%) comes from the high energy region a; > 1 GeV^. This means that the first 
nonvanishing term in the asymptotic expansion of < VV — A A >, call it c/p®, plays an impor¬ 
tant role. In the minimal saturation approximation c is predicted to be c = F^HyMl. In QCD, 
c is not exactly a constant because of the anomalous dimension of the dimension six operator 
involved. Ignoring this fact, and using the vacuum saturation approximation^ to estimate the 


value of the dimension six condensate||28||, one indeed finds a value consistent with that above, 
within a factor of two. This suggests an algorithm for systematically improving the calculation. 
Adding more resonance poles to the rational parametrization (^I]) , one can constrain the extra 
parameters so as to improve the interpolation function at both ends. On the one hand one can 
ask that more terms of the asymptotic expansion in 1/p^ be reproduced (in practice a better 
determination of the dimension six and dimension eight condensates than available at present 
is needed) and on the other hand, in the Minkowski region, one can improve the agreement 
with the experimental resonance parameters such as Ma and Fa- 

3.4 Resonance saturation of < VAP > 

Let us now investigate the approximation of resonance saturation for the three-point func¬ 
tion iy^^(p, q) (0)- As explained above, this approximation consists in constructing a rational 

® This approximation is exact in the limit Nc —> oo. In practice, however, it seems to be valid within a factor 
2-3 
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interpolation function, the parameters of which will be constrained from both the asymptotic 
euclidian region and from the physical resonance region. In the asymptotic region one should 
at least impose the constraints which are necessary to insure hniteness of the result for the 
low-energy constants Kn and K 12 . As before, the minimal version of resonance saturation has 
as resonance content the p and the Oi. Whether this is sufficient in the present situation is by 
no means obvious and must be checked explicitly. We note hrst that satishes two Ward 
identities: 




2^ {(lu (p + q)i 




(p + qf 


q^W,,{p,q) = {-BoF, 


2^ (P + g)M 

(p + g)2 ’ 


which imply that IWj, must have the following tensor structure. 


Wf,u{p,q) = -BqF^ 


(Pm + 

g 2/2 /2 


^ + F(p^g^/2)P,. + G(p^g^/2)g 


{lU 


(43) 


(44) 


with /2 = {p + qy and where P^u and are the two independent tensors which vanish under 
contraction with p^ as well as with q^, 

Ptiu = Pvq^ - P-q Q^lu = p'^qt.qu + q'^Pf.Pu - P-qPi^q,^ - P^q^Qi^v , (45) 


F and G are analytic functions of the three variables and = [p + qY with appropriate 

poles and cuts. In the simple approximate representation used here, F and G are meromorphic 
functions. 


Consider now the asymptotic euclidian region. The OPE provides an expansion valid in the 
regime where all three momenta squared p^, q^ and P are large and negative. In other terms, 
scaling p ^ Xp and q ^ Xq one obtains an expansion in inverse powers of the scale parameter 
A. The leading part in this expansion is controlled by the dimension three qq condensate and 
scales as 1/A^ (this scaling behaviour is exact as anomalous dimensions cancel out). We have 
also worked out the subleading part in this expansion, scaling as 1/A^, which is proportional 
to the so-called mixed condensate. 


lim Wf,^{Xp,Xq) 

X^OQ 

< ua^i^G^^u > 

6 A^ 


<uu> j {p^, + 2qf,)q^ 


A 2 


piqi 


+ 


( g 2/2 

p2 _|_ ^2 p2 _|_ ^2 

pH'^ qH^ 


p"^ -q^ - P 
P 2p2g2/2 


p 


flip 


pH^ 



There are also 0{as) corrections which we have not evaluated. With the minimal resonance 
content, it proves perfectly possible to match the leading part of this asymptotic expansion. 
This may be viewed as a surprise as this matching is much more constraining than that for 
the two-point function < VV — AA >. We have now three independent variables instead of 
just one and two independent amplitudes F and G. In fact, this asymptotic constraint entirely 
determines the rational approximant up to two constants, a and b 

2 ,24 ^ (p^ -q^ -P + 2a) 2 2 / 2 n _ -q^ + b _ 

^ 2 (p 2 -M 2 )(g 2 _Ml)/ 2 ’ - M^){q^ - M\)qH^ ' 



F(p 2 


Note that the effective lagrangian model used in refs. |n 
constraint. Indeed, using this model one hnds 


is inconsistent with this asymptotic 


F{p\q\P) 


1 (F^-2FvGv 

F2 (p2 - M2 )/2 


(g 2 



G{p\q\P) 


1 —2FyGy 

T! (p2 - M2 )g2/2 • 


(48) 
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Comparing with the QCD prediction in the asymptotic region, eq.(|^) one observes that the 
scaling behaviour is correct but one does not match the individual terms exactly. Worse than 
that, this model for < VAP > does not satisfy the weaker constraint to reproduce the correct 
QED divergence in the calculation of 2Kii±Ki2 (in other terms, after taking due account of the 
counterterms, the result for Kn and K 12 is inhnite in this model). The reason for these problems 
is that no vrpai coupling has been introduced. Allowing for such couplings one can recover the 
result of eg. (|T7D . If one now attempts to match not only the leading asymptotic terms but also 
the subleading ones in the expansion (^61) one hnds that this is no longer possible unless the 
resonance content is enlarged. Obviouly, for instance, in order to produce terms proportional 
to 1//^ we need not only the vr but also the tt' resonance, coupling to the pseudoscalar current. 

Returning to the minimal model, it turns out that there are further chiral symmetry con¬ 
straints as well as asymptotic constraints. These constraints concern the pion to vacuum matrix 
elements that one can deduce as residues of the pion poles in < VAP >. The minimal model 
will prove capable of obeying all these constraints and the two constants a and b will be deter¬ 
mined. Consider chiral constraints hrst: the pion to vacuum matrix element < 0|RA|7r > must 
satisfy a soft pion theorem]^ 


-Co 


(49) 


Using the rational parametrization for < VAP >, eg. (^71) , and the corresponding one for 
< VV — AA >, eg. (PD, one hnds that the soft pion theorem is obeyed provided the hrst 
Weinberg sum rule holds and the parameters a and b satisfy 


a-b= - MjFp = -(Mj + Mj) . (50) 

-Co 

where the second equality follows from imposing the second Weinberg sum rule. A second 
soft pion theorem associated with the pion to vacuum matrix element < 0|UP|7r > is satished 
without bringing new constraints. 

We are thus left with a single arbitrary constant, b. The physical meaning of this constant 
is made clear by identifying the vector form-factor of the pion. For that purpose, let us consider 
the residue of the pion pole (p — = 0 in 


lim (p - 1) = iPoip^ - 2l^){p^ - Q Fy(p^) . 

(p-iy^o 


(51) 


Here, the function Fv (p^) is the vector form factor of the pion dehned in a standard way 


<7r“|uv'= >=*r>::+p^)Pv(p^). 


(52) 


Using our model for W, we obtain the following expression for Fv{p^) 


Fvip^) 


bp‘^ 

2Ml(p2 - m2 ) • 


(53) 


It is thus tempting to determine b in order to reproduce the standard VMD form of the vector 
form factor of the pion. Amusingly, this property need not be imposed by hand here, but can 
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be deduced from the operator product expansion}]. Consider, indeed, the pion matrix element 
of the product of the vector and pseudoscalar currents, 

W^{p,l) = I < 0\TV^{x)P\0)\7r{l) > . (54) 

The OPE implies that for large euclidian values of p^, 1) behaves as 


W,(p,l)^BoFo\^ + 0{\) 


(56) 


This is to be compared with the result obtained from the rational parametrization of < VAP >: 


W,{p,l)=BoFo 


hi 

p2 


p ■ I, b 

i + hr(-2 + 


p^ ' M\' 


+ 


hi 

p2 


2 - 


Ml 


+ 0 (^) 


Clearly, this will match with the OPE result provided one takes 

b = 2Ml , 


(56) 


(57) 


which is the same value that also insures that the pion form factor satishes VMD exactly. 
Furthermore, one can work out the asymptotic expansion of the pion to vacuum matrix element 
1) when p^ -oo. 






pZ 


V QlivP ' I Puh P 

p2 \ 


p-l\ ^ p^hh -p-lp^ly 


pz 


+ 


p‘i 


+ <^( 3 ) 


(58) 


It can be checked that this result is exactly reproduced in our model provided eq.(bOD holds. 
Inserting this asymptotic expansion into eq.(3D) one verihes explicitly that the UV divergence 
cancels out with the counterterm contribution. 


At this point, we have verihed that the rational parametrization of < VAP > in terms 
of two resonance poles (and the pion pole) is capable of matching the leading terms in the 
asymptotic expansion of < VAP >, as well as the leading asymptotic terms of the two related 
pion to vacuum matrix elements, while obeying all the chiral symmetry constraints. Using 
this parametrization (see (^4}) , (^7]), (^), (^)) it is not difficult to compute the momentum 
integrals in (^^ , (|39|) and match the result with that of the chiral expansion for small values 
of external momenta. One hnds that the QED inhnities do cancel out exactly, and one obtains 
the following hnite results for the low energy constants Kn and Ki 2 '. 


^n(h) = 


1 1 


Ibvr^ 8 




+3'°sd2(z-l)3+d + 


■V 

3z{z-3)' 

JZZTf 


(59) 


where 


z = 


Ma 

M2 


^ This actually holds only for the minimal rational approximation. 


(60) 
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and 


= 


1 Ij., 3 


leTTMT'" 


^log 


V 


M 2 


^ 73 , z +1 3z 
« + 4 + 2'°s"(7rT)5- — 


(61) 


The dependence npon the chiral renormalization scale fi agrees with that derived in ref. 

(see also |@) in the gange ^ = 1. As was anticipated, one observes that Ku and K 12 also 
depends on the QCD renormalization scale /iq. The resnlt on this dependence agrees with 


ref.|^ concerning it'll bnt not it' 12 . We stress that once the assnmption of satnration from a 
single mnltiplet of vector and of axial-vector resonances is made, the result for Ku and it '12 is 
uniquely determined, in much the same way as it was for the 0{p^) constant C, in terms of 
just Ma and My. We also note that once the expression for K12 is obtained, the three-point 
function < VAP > may be used for getting K 13 via (|35D . We checked that the same result for 
Ki 3 which was found before from a two-point function (i.e. eq. (^7|) , plus resonance saturation) 
was indeed recovered. 


3.5 Phenonenological implications 

Our construction of the Green’s function < VAP > has phenomenological implications 
which are interesting to consider in order to assess the reliability of the result for Ku and K 12 . 
Let us hrstly consider the predictions in the region of low momenta and then, further away, in the 
resonance region. For small momenta, comparing with the chiral expansion of < VAP > (which 
can be found in ref.]^), one hnds that the model reproduces the vector meson dominance 
formulas for the constants ll{My) and ll{My) which are known to be reasonably accurate 
Obviously, however, the model does not generate the logarithmic singularity caused by the pion 
loop. It is perfectly feasible to improve on this by taking a more sophisticated form for the 
vector meson propagator. However, this will complicate the computation of the photon loop 
integral while bringing corrections which are subleading in the large W counting, and should 
thus be rather small. Concerning the resonance sector, for the p meson, the model is found to 
embody a reasonable value for the tttt decay width. For the a\ resonance, the model predicts a 
value for the pvr decay width of the order of 200 MeV, which is within 50% of the experimental 
value. A surprising result emerges for the radiative width, ai —>■ yTr. Extracting the decay 
amplitude from < VAP > using the reduction formula, one obtains the following expression 
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T(ai(g) ^7(p)7r) = {e.pe'.q - p.qe.e') 


MaFa 


2M^ - Ml 

2Ml 


(62) 


Because of the factor 2My—M|, this amplitude is strongly suppressed and one obtains a value of 
a few tens of KeV for the width, which is one order of magnitude smaller than the experimental 
value obtained in ref. [^. Clearly, since the Oi pole is rather far from the integration region 


relevant for the evaluation of the EEC’s Ku and K 12 it is plausible that these constants could 
be reasonably well evaluated using the rational approximant for < VAP > and at the same 
time, the properties of the oi need not be too precisely reproduced by the approximant. Still, 
a mismatch by one order of magnitude would be disturbing. A closer look at the literature, 
however, reveals that a small value for the radiative width is not ruled out. Indeed, a recent 
photoproduction experiment has found no trace of Oi production ( while observing very 
clean evidence for 02(1320) production). In this experiment, the photon is on shell, which is not 
necessarily exactly satished in the experiment of ref. which is a Primakov-type experiment. 
Our amplitude is very strongly energy dependent and it is nearly vanishing only for exactly 
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massless photons. At any rate, if the Oi radiative width turns out to be unsuppressed, this 
would be clear evidence for the necessity of including further resonances in the parametrization 
of < VAP > and would invalidate the above estimate of Kn and K 12 . On the contrary, 
conhrmation of the suppression would be a nice experimental support for these estimates. 

3.6 Some remarks on Ki to Kq 

Let us discuss here some aspects of the convolution representations for the constants Ki to 
Kq, with some emphasis on the combination + Kq which is of interest for DT violation. In 
order to illustrate the procedure, consider the correlator 


Traftcd 


ip) 


= I < 0|T [Al{x)A»^iy) - V:ix)V^iy)) Q'yiz)QU 0)\0 > (63) 


We will concentrate in the following on the hrst tensor component, LIq^'^'^. Using the chiral 
lagrangian, we can compute its chiral expansion, which involves a photon loop as well as a pion 
loop. Choosing flavour indices in order to isolate + Kq, one obtains. 


Klip) + Klip) - -KUp) + 3KUp) 


(64) 


+ 


1 3 


167r^ 4 


log 




p 


^ + - 
2 


log 


, eioge 

h' e-i 


+ 




-log 


—p 

p^ 


+ 1 


+ Oip^ 


This expression allows one to deduce the scale dependence of K^ + Kq, utilizing the previous 
results for K 12 (eq.(pT|)) and iCia (eq.(P7|)), 




3C' 


( 66 ) 


which reproduces the result of ref.[IJ] with the new feature that it is found not to depend on 
the gauge parameter ^ (we note in passing that that the correlator < VaVyQvQv > is hnite at 
one loop, which immediately implies that d/dpiK^^ip) + 2Kl^ip)) = 0 ). In a similar way, it 
is easy to pick up two sets of flavour indices such that the chiral expansion of the correlation 
function will involve the combinations Ki + K 2 , 


3C 

W 


nr“(p") = |AT(p) + K'M + -(A';(p) + kUi^)) + 

and —2K^ + K^, 

{o Qn 

-2Klip) + Klip) + ^iKlip) + Klip)) + ^ 


log- 


—p 

p^ 


+ 1 


( 66 ) 


log^ + 1 

K 


(67) 


The next step is to express the correlator above (^) by taking functional derivatives 
QED and QCD actions. One observes, hrstly, that the counterterms make a contribution 
has the following form. 


of the 
which 


KTiP)\counter = *^2(n“f(p/2) - n“^(p/2)) . (68) 
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From this, one deduces the finite convolution representation 


Tt1425 


(p) 


= —I 


d^k 


D^,y{k) 


n' 


-iF^-Zo 


9aP 


PaPp 

p 2 


Ai!44y2y5(p, k) 

(l + 0(p^)), 


n 


aPfiu 

yly4y2y5 


(P. k) 


(69) 


in terms of the QCD 4-point functions, 


yiaPliu 


{p,k) 


I < 0|TA“ (x)p4j(|/))V;"(^) 1/^(0) |0 > 


(70) 


and the similar correlator with y4^ replaced by 17“, V^. 


The analogous representations for and have no contribution from the coun¬ 

terterms implying that the convolution integrals must be finite. One must extract the tensor 
component from the four dimensional integral and one obtains + Kq using eq.(M). 


The part proportional to — 1 in the convolution integral above (|^ can be worked out 
in a model independent way. Indeed, the contraction of k^k'^ with the 4-point function can be 
expressed in terms of the 2-point function < VV — AA > using the chiral Ward identities. In 
this manner, the dependence upon ^ of ng^^®(0) can be exhibited as. 


nr(0) = nr(0)k.„ + (gjTj jiog^ + log^ - 4 + log^ |pv(n - pVx)|| 

(71) 

in which the expression of the counterterm constant Z 2 (see (pTf) ) has been used. This formula, 
together with the expressions of K 12 and Ki^ (0). o allows one to show that the combination 
+ Kq is independent of the gauge parameter Furthermore, the dependence on the QCD 
renormalization scale /xq is found to be entirely absorbed by that of the parameter K 12 such that 
Kq + Kq is independent of /iq. In a similar way, one can also show that the two combinations 
K 1 + K 2 and —2KQ+Ki are gauge independent and po independent. Correlation representations 
for all of Ki to Kq can be obtained by further considering the correlators < A'^A^^Q'iQj^ > 
from which one can isolate K 2 , K 4 and Kq. The derivation exactly parallels the discussion 
of the correlator (^) and will not be detailed. Let us simply mention that there will be no 
contribution from the counterterms, implying that the convolution integrals must be finite. 
Again, the dependence on the gauge parameter ^ can be worked out explicitly and it turns out 
that the six parameters Ki to Kq are independent of the value of 


In order now to obtain phenomenological estimates, for instance for Kq + Kq, we must 
construct a resonance saturated model for the QCD 4-point function < [AA — VV)VV >. 
The asymptotic constraint to be satisfied in this case is that the convolution integral must be 
finite in the particular gauge = 0 (Landau gauge). The model must also obey a set of chiral 
Ward identities. The simplest way to implement these in the present situation is to first write 
down a chirally invariant effective lagrangian for the resonances. For this purpose, one may 
use, for instance, the convenient method advocated in ref.f^. However, it will be necessary to 


investigate more couplings than were considered in this reference (where only those couplings 
actually contributing to the LEC’s Li were discussed). We already pointed out that in order 
to fulfill the asymptotic constraints associated with < VAP >, vertices associated with aipvr 
couplings were needed. Here, we can also expect contributions associated with ujpn coupling 
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(its contribution to AM|- was considerered long ago 0) but possibly also from hiojn coupling. 
Also four-particle couplings and exchanges of scalar or tensor mesons could be considered. It 
is clear that a rather extensive study is necessary in order to carefully investigate all these 
possibilities. This study will not be undertaken here. We have simply performed the exercise 
to reevaluate the LEG iGis, which can be isolated from the 4-point function < v^v^QlQi >, 
using the resonance effective lagrangian of Ecker et al.[^. The number of graphs to consider 
is strongly reduced by working in the Landau gauge (the gauge dependence can be worked out 
explicitly as above). The calculation is much more tedious than the one based on a two-point 
function in sec.3.1, but one indeed recovers exactly the same result. 

4. Flavour symmetry breaking and the constants 

4.1 The constants Ky and Kg 

We can isolate the contributions of the two constants K-j and Kg by taking the derivative 
with respect to of the Q-correlators < Qv^Qv^ > and < Qa^Qa^ > already considered. 
After a simple calculation, one obtains from the chiral lagrangian, 

= 2BoFS |2(AT(^) + K'M) - + log^)| + 0(m,) (72) 

and 

Wn« (0) = 2 B,FS i^2(K'{p) - K’M) + + log | + 0(m,) . (73) 


These expressions display the scale dependence of Kj and K^. In particular is hnite and 
scale independent. One may then use the convolution representation for the Q-correlators Ilys 
and n ^3 in terms of QCD correlators but this does not lead to an expression which can easily 
be evaluated. However, it is not difficult to show that the derivatives with repect to the strange 
quark mass of the isospin one vector and axial-vector spectral functions are suppressed in the 
large A), counting. This is because the graphs which contribute involve one extra (strange) 
quark loop with respect to the leading order graphs. The fact that Kj and Kg are subleading 
in Nc could, of course, have been suspected from the double trace structure of the corresponding 
lagrangian terms (see ref.for more details on this aspect). One may assume, then, that Kj 
and iGg are suppressed for scales /r ~ My and infer the following estimate. 


^7-0, ir8^(/i)^ 


c 


167r2 2F^ 


log 


K 

/i2 


(74) 


4.2 Constraints on symmetry breaking spectral function differences 

In order to discuss the constants Kg and Kiq, we consider the correlator 

An^(/) = ^ - f < 0\T{QUx)QU0) - Qv(a^)Qv(0))|0 >un • (75) 

2(m — nis) < uu> J 

and a similar correlator AH® with V replaced everywhere by A in the above formula. The 
normalization factor, in front of the integral is introduced for convenience. The subscript lin 
means that we must keep only terms linear (or logarithmic) in the quark masses and drop the 
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quadratic or higher order terms (the reason for this will be given below). Using the dehnition of 
the Q-currents in terms of the QCD vector or axial currents and the spectral representation of 
the correlator Yiysi^q) (see (^)) and the analogous definition for the following spectral 

function difference will appear, 


Apy(a;) 


-3 

2{m — nis) < uu> 




(76) 


Again, we can also define Apa{x) by replacing pya by everywhere in (|7BD. In the latter 
case, as we are away from the chiral limit, there is a second spectral function a a (see (|^) ) to 
be considered and we define 


AaA(x) 


-3 

2(m — m<j) < uu > 




(77) 


With these definitions Apy and ApA must satisfy sum rules analogous to the two Weinberg 
sum rules. Indeed, at large —q^ , the leading term in the asymptotic behaviour of the QCD 
correlation function difference < — U®U® > is given by| 


lim 





< uu > 



(78) 


This implies the two sum rules, 

POO POO 

/ dx Apy{x) =0, ly = — dxx Apy{x) = 1 . 
JO n Jo 


(79) 


For these sum rules to hold, it is essential to drop the quadratic mass terms in the definition of 
Apy{x). Otherwise, the first sum rule would still be valid but the second one would diverge. 
The first of the above sum rules was established long ago|l3l| (and is customarily referred to 
as the DM0 sum rule) while, curiously, we could not find a trace of the second one in the 
literature. Analogous sum rules also hold for the axial current spectral functions 


POO POO 

/ dx Apa{x) = 0, Ia = — dxx Apa{x) = 1 , 
Jo IT Jo 


(80) 


together with 


1 

TlJO 


POO 

- / dxAaA{x) = 2 . (81) 

'Jo 

The latter sum rule is in fact saturated by the pseudo-Goldstone boson contributions to AaA, 


AaA^x) = 2 


M^S{x - Ml) - Ml6{x - Ml) 


M2 - M2 


+ Atr^)!-) 


(82) 


Because the divergence of the axial current is linear in the quark masses, the piece AaA{x) is 
of higher order in the quark expansion than the first term in (|5^) and must be dropped. 

Instead of the neutral vector currents and U® (or the corresponding axial vector currents) 
one could equally well employ charged currents and 


jrud 






(83) 
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(and their hermitian conjugates). In particular, the sum rules (|79D and (^0]) hold unchanged 
for the spectral function difference of the charged vector current Apy defined as follows, 


^Pv{x) = 


[rh — m.) < uu > 


(84) 


and the corresponding definition for the axial-vector case. The charged vector currents are not 
conserved but this effect is proportional to the quark masses and can be ignored here. 

4.3 Sum rules for Kg and Kiq 

Let us now return to the Q-current correlator Any(p^) dehned in (|75|). As before, we hrst 
compute this correlator from the chiral lagrangian, up to chiral order four. Setting = 0, one 
obtains: 


An^(O) = 8 


where 


Z°(p) = 




1 3 


Z°(p) 


-I- 0{m, 


IGtt^ 4 

C Mi log(Mi/p2) - m 2 log(M2/p2 


<?/ ’ 


^ 0 ^ 


Mi - Mi 


( 86 ) 


( 86 ) 


is the pion (and kaon) tadpole contribution. In terms of QCD currents, we obtain the convo¬ 
lution representation 


An«(o) = 


-3 


d^k 


- ni"s(fc)) {-iD^,{k)) + Zs-Zg. (87) 


2(m — rris) < uu > J (27r)'^ 

Next, inserting the spectral function representation for the correlator < > yields 

An^(o) = 


2 roo p 2 

dxx/\pv{x) V{e)+ \ogA'K + --\ogx 
lOTT^ TTJO 3 


-]r Zg — Z‘ 


2 • 


( 88 ) 


Some remarks are in order here concerning the ultraviolet divergence. From the large momen¬ 
tum expansion of the correlator < >, one observes that there are two terms 


which will lead to a divergence in the photon loop integral eq. (|87|) . 


- ny8(p) = 



^ (m — ms) < uu > 


1 + AA + ,,, 

37r 


+ 0 (- 


pD 


(89) 


(retaining only terms linear in the quark masses). The first term in the square bracket in eq.( 
generates a 1/e pole in dimensional regularization. This pole appears explicitly in the spectral 
representation (^8]). Using the expression of the QED counterterms Zg and Z 2 dO) one easily 
verihes that the 1/e inhnities cancel exactly. This cancellation requires that the second sum rule 


([79|) be satisfied. It is important to keep this point in mind in phenomenological applications 
based on resonance saturation. In particular, it is not consistent to ignore flavour symmetry 
breaking in the resonance sector as has been done in ref. . If one does not impose that this 
sum rule be obeyed, the result will be inhnite. Now the second term in the square bracket (|89D , 
proportional to ag, also generates a divergence in the photon loop integral, of the form loge 
in dimensional regularization (such a smooth singularity is the result of renormalization-group 
improvement). This divergence should be removed by minimal subtraction. This means that the 
integral J^dxx logxApv{x) appearing in eq.(|^) is, in fact, divergent (the same holds true for 
the axial spectral function Ap^)- Since the integral J^dxxApv(x) is convergent, this suggests 
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that the spectral function difference behaves as Apy ~ log^x asymptotically. This remark 
can be used to minimally subtract the loge singularity in the spectral representation, which 
should be done in practice if one were to use realistic spectral functions. The simple models 
which will be considered below do not lead to such singularities, so we will ignore this subtlety 
in the following. 


In a similar way, one obtains for the corresponding axial current correlator, firstly from the 
chiral lagrangian 


An^(O) = 8 




e M 2 log(M2/p2) _ m 2 log(M2/p2) 


167r2 


4(M? - Mi) 




167r2 4 


and secondly from the QCD action 


AnS(o) = 


-3 


d^k 


2 (fh — m,s) < uu> J (27r)‘^ 




+0{mq) 

(90) 

2 ( 91 ) 


2 roc p 

— / dxxAp^ix) r(e) + log47rH-logo: 

3 


167 r 2 TT Jo 

I I roo 


-J dxAaA{x) (^ + 3)(r(e) + log 47 r + 1 - logx) 


lOvr^ TT Jo 

Again here, the sum rules (|80|) and (^ID for ApA and Acr^ ensure the correct cancellation of 
infinities. One notices that the pion contribution in the photon loop in eq.fpOD appears also 
in eq.(^) in spectral representation so that this contribution cancels out in the expression for 
Kq — Kiq. One also expects that the chiral logarithms in the function Z^{p) in eqs.(^) and (|90|) 
should cancel out with similar terms present in the spectral functions but we will not attempt 
to elucidate exactly how this happens and we will keep this contribution (which is subleading 
in the large counting) as it is for the moment. 

It is now easy to deduce the sum rules for Kg and Kiq. Defining the two integrals 

poo OP poo 

Zv= dxxlog—jApvix), Za= dxx \og—^ApAix) (92) 

JO My JO Mv 


'y 


where ApA is the part of the spectral function with no pion contribution. 


Apa{x) = 2 


F^J{x - Ml) - F^^5{x - M2) 
(M2 - M2)F2 


+ Apa{x) 


(93) 


(and the worrisome asymptotic tail related to is assumed to be removed). We hnd the 
following result for Kg\ 


Klip) = 


1 1 


p 


167r2 8 


-U\og^-i--{Zv + ZA) 


ho 


and for A'lo, 

^[o(h) = 


1 1 


167r2 8 


+ (^ + 3) log + ^ + — Zy) + 1 + QZ^{p) 


(94) 


(95) 


The same expression for the constants Kg and Kig holds in which the spectral function dif¬ 
ferences related to the neutral currents Apv and ApA are replaced by the charged ones Apy 
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and Ap^ (see eq.(F^)- The only modifieation to eqs. (p^ and (|9^ ) is that the tadpole funetion 
Z^{p) must be replaced by Z^{p), with 


Z+(p) = 


c 




Mi) 


3M2 log + 2 Ml log - 5M^ log 


M2 


Mi 






Ml 

/i 2 


(96) 


4.4 Difficulties of minimal resonance saturation 


As one can see from eqs. (|9^ , (|95D and from (|108|) , (|109|) below, in order to make a definite 
prediction concerning the violation of Dashen’s theorem one must be able to evaluate, in a 
reliable way, the the integrals Za and Zy defined in (|9^. A priori, we will envisage to estimate 
these integrals in the approximation of resonance saturation, retaining the contributions of the 
lowest-lying vector meson and axial-vector meson octets. One must first verify that the spectral 
functions /S.pv{x) and Apa{x) approximated in this way obey the sum rules (^) and (|80[) with 
acceptable values of the resonance parameters. The first sum rule for Apy and the equivalent 
sum rule for Apy imply, respectively. 




Fp+ — Fk* 


(97) 


where the couplings of the neutral vector mesons are defined with respect to the electromagnetic 
current, i.e., < >= MyoFyoe^. Assuming exact isospin conservation allows one to 

extract the matrix elements of the currents and V^. Using the e+e” branching ratios, one 
obtains 

Fpo = 152.9 ±4, F^ = 45.9 ± 1.0, F^ = 79.1 ± 1.5 (MeV) . 

If one now extracts the value of FpO using the first relation (^), one obtains FpO ~ 158. 


(98) 
This is in 

rather reasonable agreement with the experimental value. We can also extract the experimental 
values of the charged decay constants from r decay data. Using the most recent compilation 
one gets 
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Fp+ = 146 ± 2, 




= 151.5 ±4.7 


(99) 


The values of Fp+ and Fk*+ turn out indeed to be nearly equal as is demanded by the DM0 sum 
rule for Apy saturated by a single resonance multiplet. This is not a completely trivial result if 
one thinks that typical coupling constants such as F^^ and Fk differ by 20%. In conclusion, we 
seem to find that minimal resonance saturation is a good approximation as far as the first sum 
rule for Apy is concerned. Let us now examine the analogous sum rule for the axial currents. 
Experimentally, the only accessible data concern the charged currents. Using the axial version 
of the DM0 sum rule for together with resonance saturation, implies 


^Ki{1270) ~F 


^ al ^77 


Fi 


( 100 ) 


In this relation and Fk are, of course, rather accurately known. Concerning Fa^, an estimate 


can be made under the hypothesis that the decay amplitude r —> 7r“7r°7r° (the most recent 
PDG[^| value for the branching fraction is 9.27 ± 0.14%) is dominated by the Oi resonance 
and proceeds via Oi —> pn (in other terms, we take r(r —>• a^i^r) = 2r(r —>• 7r“7r°7r°z/.r)), which 
gives Faj = 165 ± 13 MeV with M^^ = 1230 ± 40 MeV). The first published data for the tan 


decay rate into iLi(1270) and iLi(1400) by the TPC /27 collaboration [P5| have rather large error 
bars. Much more precise results obtained by the ALEPH collaboration were presented very 
recent lyp 6 


Transcribed into coupling constants these results read 


Fi 


iCi(1270) 


= 1232 ± 63" 


p 2 

-^if2(1400) 


= 972 ± 682 (MeV^) 


( 101 ) 
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The errors are certainly larger than in the vector case, still it appears plausible that the axial 
DM0 sum sule is reasonably well saturated from the lowest lying resonance contributions. 

Let us now turn to the second sum rule that /S.py and must satisfy ([75D, (|5D|). Consider 
the vector currents hrst. Computing the integral ly using minimal resonance saturation and 
enforcing the hrst sum rule gives 


- M}) 

F'^yp, - Ml) ■ 


( 102 ) 


The second sum rule states that ly = 1. Unfortunately, this relation completely fails to be 
satished as, using experimental values, one hnds ly ~ 5.4. Neither experimental uncertainties 
nor the use of the narrow width approximation can explain such a large discrepancy. One source 
of uncertainty stems from the requirement of expanding the numerator to linear order in the 
quark masses. We have assumed that the term linear in the quark masses dominates the chiral 
expansion of the vector meson masses, as is suggested by the success of the Gell-Mann-Okubo 
mass formula for the vector nonet. The only possible explanation, then, is that the integral 
ly picks up signihcant contributions from the energy region above 1 GeV, implying that the 
simplest resonance approximation method appears to fail in this case. As has been discussed 
in sec.3 this approximation consists in performing a rational interpolation of the correlator 
< > with two poles. In the present case, the interpolation is not capable of 

correctly matching both the region of large and the region of small In order to improve 
it, we must increase the number of poles. If one nevertheless insists on employing a single 
resonance multiplet, then one must use an unphysical \a\\ie for the product F^{Ml-Ml) such 
that ly in ( |102|) satishes ly = 1 { recall that if the sum rule were not satished the result for 
Ag + Aio would be inhnite). Doing this we hnd the following value for the integral Zy (|104|) 
which occurs in the expressions of Ag and Aig: 


= Iy = l . (103) 

Equality of and ly results from using the narrow width approximation, together with the 
hrst sum rule, and expanding linearly in the quark masses. This result, however, is not stable 
against inclusion of higher mass resonances. Indeed, let us include one additional multiplet of 
vector resonances. One can ht the value of i{Mp, — M|/) as well as the splitting F^, — A|/ 
in order that the hrst two sum rules be obeyed^ while keeping physical values for Fy and 
Mp — M^. This certainly improves the validity of the interpolation function for small values of 
the p^. Whether the values of the new couplings Fp, F^i are at all realistic cannot be decided 
since not much is experimentally known about these quantities. This improved model yields 
the following evaluation of Zy: 

®This is possible provided — 3F^ — 3F^ is negative which, fortunately, happens to be the case. 
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(104) 


Zr^pr ^ ^ ^ 


3(F^^ - - 3F|)M^ 

2(M2 - Ml)F^ 


'M^, , m 2 


M 2 


^1 ^V' 


1 - 


(M 2 


Mi)F 2 ' 


m - MI)F^ 


~ -3.9 


taking My = 1.6 GeV, which differs considerably from the estimate based on one resonance 
multiplet. 


Turning to the axial-vector sector now, let us evaluate the integral I a in (^OD under the 
approximation of minimal resonance saturation. Assuming the first sum rule to be satisfied, 
one obtains. 


Ia = 


(Mi 


Ml)Fi 


+ 2M'i(Fi - Fi) 


+ 2 , 


(105) 


and IA should be equal to one according to the sum rule (^). Using experimental numbers, 
one finds instead J/i ~ 6 ± 2. If one uses for Ma and Fa the values from the Weinberg sum 
rules (i.e. Fa = 122 MeV, Ma = 966 MeV with our choice of Fv and My and ignoring the 
difference between Fq and F^^) rather than the experimental ones, then one obtains I a — 4. 
The approximation of minimal saturation is again found to be in conflict with the second sum 
rule, if one takes physically reasonable values for the resonance parameters. 

At this point, the conclusion would be that it is not possible to obtain a reliable estimate 
of both Zy and Za (and consequently of the two LEC’s Kg and A'lo) at present, because the 
spectral integrals are too slowly converging. However, if one is only interested in the amount of 
Dashen’s theorem violation one really needs to be able to estimate A'lo, because Kg happens to 
be multiplied by M^ and will contribute very little. Now the convolution representation for Kig 
involves the double difference Hys — Hy^s — Hys -|- H^s. The corresponding spectral functions 
must satisfy two sum rules obtained by combining ([79|) and (^OD 


rOQ poo 

/ dx{Apy{x) — Apa{x)) = 0 , / dx x{Apy{x) — Apa{x)) = 0 

Jo Jo 


(106) 


The second sum rule implies that convergence of the spectral representation should be faster 
than it is for Hys — Hys and H^a — H^s individually. In fact, one can observe that the second 
sum rule in ( p.06|) , which can be stated as I a = Ivi is approximately satisfied in the minimal 
resonance approximation as one can check from the values of I a and Iv obtained above. Ex¬ 
perience with the Weinberg sum rules and the corresponding DGMLY sum rule for C suggests 
that minimal resonance saturation is likely to be a reasonable approximation as far as Kig is 
concerned. In this approximation, the value of the difference Za — Zy to be used in eq.(^) for 
Kiq is given by 


Za — Zy = 


F^AMl - Ml) 


Ml 


2 Ml{Fl - Fl) + F2 (m; - M2 ) log ^ 


- 2 1 -h log 


Mi, 


(107) 


where the sum rule I a — /v = 0 has been used to simplify the expression. 
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5. Consequences for Dashen’s theorem violation 

Up to corrections which are quadratic in the light quark masses or quartic in the electric 
charge, the contribution to the mass difference proportional to can be written as 

a sum of six terms |^, as follows: 


AM|- = M^+ — M^o |em — to + + ^2 + ^3 + ^4 + ^5 


(108) 


with 


tn = e 


U = 


,2C 


^ 0 ^ 


IOtt^ 


3M|- log 


Mi 


4:Mf,+ 


Ml 


u = 


.lec r 




{Ml + 2Ml)Ll + Mia 


h = [4A'J + -(K; + A-) 

U = eHll [SKI + 8(AT„ + A-;,)] 
t5 =-e^MjhA-+ A-J) . 


(109) 


M2log^ + 2Milog 


MM 


One notices first, using expressions (p4|), (^Sj) for Kg and Kig and (^9]) for Ku and the discussion 
in sec. 3.6 on that the dependence on the QED gauge parameter ^ drops out, as it 

should. In fact, the sum to + ti + ^2 + ts + and are separately gauge independent. The 
dependence upon the QCD renormalization scale fig drops out in the term but not in the 
term tg. This indicates that the electromagnetic contribution to the mass difference 

is not, strictly speaking, an observable quantity: only the full mass difference is independent 
of /To- Fortunately, this term tg makes a negligibly small contribution to the mass difference, 
being proportional to Ml. In order to estimate numerically the value of AMl one needs, 
essentially, an estimate of the three constants Kig, Ku and K^ + Kg ignoring the terms which 
are proportional to Ml or suppressed by the Zweig rule. The estimate of Kig based on resonance 
saturation is contained in the expressions (^) and (|107|) and a similar estimate for Ku is 


given in (M). We were not able to evaluate Kg + Kg but one notices that this combination 
appears in the expression for AMl with a much smaller numerical coefficient than Kig or Ku- 
Assuming that Kg and Kg have the order of magnitude typical of 0{p‘^) LEC’s the lack of 
precise evaluation generates a rather small uncertainty. Based on these estimates one obtains 
the following decomposition of AMl: 


M^ 


K+ 


■M^ 


\ em 


(M2+-M^o) {l.02+ 1.13 -0.20+[-2, 2]10-2 +0.80+ [-0.3, 0.3]} . (110) 


where the six entries correspond to the six contributions displayed in eq. (|109|) . Indi¬ 

vidual contributions are scale dependent: they are shown at the CHPT scale /i = My and the 
QCD scale fig = 1 GeV. The values of Ma and Fa wherever they appear are taken from the 
Weinberg sum rules. The contribution from tg (fourth entry in eq. (|110|) ) has been estimated 
to lie in a range, assuming that Zy could take values between ±4 (see sec.4.4). The last entry 
corresponds to the contribution of Kg and Kg assuming that each of these LEC’s lies in the 
range ±l/167r^. 
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It is interesting to verify the stability of the resnlt if one drops systematically all terms which 
are known to be snbleading in the large Nc connting. This is easily accomplished in eq.( |lU 8 |) 
knowing that Fg, and C are 0{Nc), is 0(1), Fg, Kg, F[g, are also 0(1) while 
Fg is 0(l/iVc). The variation of LK^) with respect to the scale /x is snbleading in bnt, in 
practice, it is not negligible since Lg varies by a factor of two as fi is varied between and 
Mp. This raises the qnestion of which scale shonld one choose for performing a leading large 
calcnlation. We bypassed this difficulty by expressing directly L 5 in terms of experimental data 
(i.e. Fk and Ft^) dropping the contributions which are subleading in W in the chiral expansions 
of Fk and Ft^. In this way, one obtains the scale independent estimate. 


L5 


_ p'i 

8(Mi - M2) 


~ 2.310"^ 


( 111 ) 


Using this result, and dropping all subleading terms in W we obtain the K~^—K^ mass difference 
in the following form, 

{ 1.02 + 0.75 - 0.51 + [-1, 2]10-2 + 1.20 + [-0.3, 0.3]} . (112) 


Interestingly, though individual terms are substantially modified when removing these contri¬ 
butions which are subleading in Nc, the overall result is reasonably stable. The contribution 
of + Kq still remains to be explored in detail, before a quantitative estimate of AM]^ can 
be provided by this approach. Still, the results above are strongly suggestive of a rather large 


amount of violation of Dashen’s theorem, similar in magnitude to the results of | 
or to the calculation performed on the lattice in the quenched approximation. 


H, M 


The cause for the different result obtained in ref. [|^ (negligible DT violation) can be better 
understood with the help of a recent preprint by the same authors]^ in which they have com¬ 
puted, from exactly the same resonance lagrangian, the contributions to each of the constants 
Ki. Firstly, their neglect of flavour symmetry breaking in the resonance sector amounts to set¬ 
ting KIq{Mv) = 0. Doing this bring the contribution of to AM]^ down from approximately 
one to 0.3, using our value for Ku. However, their estimate ^ of Ku is much smaller than 
ours such that ^4 is practically zero in their evaluation. As mentioned in ref. [^, a remainder 
contribution is needed in order to restore the correct chiral scale and QCD scale dependences. 
Comparison with our results shows that this remainder is not negligible, and may even be 
dominant as in the case of Ku. 


We note that neither [^, or 0 have included enough terms (such as vrpai coupling, 
see sec 3.4) in their resonance lagrangian such as to satisfy the asymptotic conditions needed 
to evaluate separately Kiq, Ku and K^ + Kq. A less constraining set of two sum rules were 
implemented in ref. which are derived by requiring UV hniteness of AM]^ (i.e. no contri¬ 
bution from QED counterterms), which holds in the limit = 0. This is obviously a 

minimal requirement for any realistic model but it is not satished by the estimate of ref. . 
In ref.j^l these two constraints are enforced by adjusting the two values of the axial-vector 
coupling constant Fk^ (assuming Fk 2 = Fk^) and the coupling constant of an excited 
axial resonance which they find necessary to include. These two sum rules are reminiscent of 


®The correspondance of their results with our formulas is as follows: one must drop the counterterm contri¬ 
butions and define the convolution integrals using the chiral MS prescription. The QCD correlators must, of 
course, be computed using their resonance lagrangian. In this way one immediately recovers that K[q = 0 and 
we also checked their results for Ku and K^. We could not, however, reproduce their result for K 13 . 
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those obtained above, eq. (|106|) , in connection with the parameter i^io, the difference is that 
they pick np additional contribntions associated with the combination Si^n — 4/3 (i^s + Kq) 
which occurs in the expression of AM]^. We recall that the two sum rules (|106|) did not seem to 
require significant contributions from excited resonances. A better experimental determination 
of the coupling constants Fk, and Fk2 would help to clarify these issues. We note also that the 
calculation of ref.|^^ includes contributions which are quadratic in rUg, it would be of interest to 
know the relative importance of these contributions. Our result for Kio is in reasonable agree¬ 


ment with that of ref.|T2|: dropping the tadpole contribution, setting /i = /xq = 0.7 GeV 
and .^ = 1, we obtain = 5.2 10“^, to be compared with the value = (4± 1.5)10“^ 

quoted in this paper. Finally, it would be desirable to estimate the error on the evaluations of 
C, Kiq, it'll and it's -|- Kq. In the case of C and Kiq one can use experimental data from r 
decay directly into the sum rules and this should allow for an estimate of the error as well. For 
the others, one must rely on resonance saturation and the reliability of the estimates can be 
assessed, in principle, by performing evaluations using higher order rational approximants (i.e. 
including additional resonance multiplets) and checking the stability of the result. 


6. Conclusions 


We have discussed an approach to the question of DT violation (and more generally, to the 
evaluation of the low-energy constants Ki which parametetrize all electromagnetic effects at 
order O(e^p^)) which is a direct generalization of the classic sum rule of Das et al.|I|]. Instead 
of the technique of current algebra, as used in []^, we have made use of the chiral lagrangian and 
the charge spurious. This technique is rather powerful, as one can judge from the simplicity 
of the rederivation of the DGMLY sum rule in sec. 3. As an immediate consequence of this 
method, one hnds that all the iF/s can be expressed as a convolution of a QGD n-point function 
(with n=2,3 and 4) with the photon propagator. We have displayed explicitly the contributions 
of the QED counterterms to these constants which allows one, if one wishes, to perform the 
calculation using an arbitrary regularization and renormalization scheme. In the MS scheme, 
the dependence of the parameters Ki to Ki^ on the gauge parameter ^ and on the QCD 
renormalization scale fio was entirely determined. 

A current prejudice, based in part on the phenomenology of the usual O(p^) constants 
Li (see [^) is that the values of the constants iF* (and, as a consequence, the amount of 
DT violation) should be essentially controlled by resonance physics at a scale of 1 GeV. The 
great simplicity and the accuracy of the DGMLY expression(^ certainly make it worthwhile 
to investigate this prejudice in more detail. This assumption concerning the role of the lowest 
lying resonances amounts to approximating in a minimal way the QGD correlation functions 
(which occur in the convolution expression of the constants Ki) with rational functions having 
the corresponding resonance poles. The parameters of these resonances are then subject to 
stringent constraints which can be expressed either in terms of sum rules that the exact QGD 
correlation function satisfies or in terms of matching conditions at asymptotic momenta. A 
well known example is the correlation function < VV — AA > in the chiral limit, which leads 
to the two Weinberg sum rules: these can be saturated to a reasonable level of accuracy by the 
p and the Oi resonances. This model correlator leads to the DGMLY formula and, as we have 
shown, to a similar expression for the constant Ki^. We have displayed a generalization of this 
construction to the QGD three-point function < VAP > in the chiral limit. Again, we have 
shown that a model containing as poles only those of the p and ai resonances (together with 
the pion pole) can exactly match all the relevant asymptotic constraints. The properties of 
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the p and the Oi resonances are then predicted to be in reasonable agreement with experiment 
provided the oi radiative width is indeed strongly snppressed as has been suggested by a recent 
experiment [0^] . These asymptotic constraints cannot be avoided as they ensure that the QED 


ultraviolet divergencies are exactly cancelled by the QED counterterms. We have shown the 
necessity of not ignoring the vrpai couplings (as was done in earlier work) if the contribution of 
the Oi resonance is to be correctly evaluated. Under the assumption of dominance of the lowest 
resonance multiplet in the vector and the axial-vector sectors, the constants Kn and K 12 are 
uniquely determined in terms of Mp and The constant K 12 does not participate in DT 

violation but it arises in other interesting isospin violating phenomena. For instance, as shown 
in ref.[0, it is the only constant Ki which appears in the Kl^ form factor and it also 

appears in the decay constants and Fk (note that these constants are no longer physically 
meaningful quantities in the presence of electromagnetism: they are gauge dependent and in¬ 
frared divergent). These estimates can in principle be improved by using rational approximants 
of higher order, i.e. including more resonances. 

A second class of EEC’s that we have considered are iCg and Kiq. We displayed sum 
rule expressions in terms of the spectral function differences pys — pv& and — p^s or, 
alternatively, pyud — pv^^ and p^ud — p^us which are quantities measurable from r decay data. 
In principle, then, ATg and Kiq could be evaluated without any recourse to resonance saturation 
models using data from e’''e“ scattering as well as r decay. A number of classic chiral sum 

in 


rules have been already explored in this manner 


At present, however, there are rather 

but the data 


large uncertainties as far as pys (or p^s) is concerned as is discussed in ref. P? 
is constantly improving. In this sector, the constant Kiq is expressed in terms of the double 
difference pyud—pyus — {pA’^d,— Pa'^‘’) for which simple resonance approximation can be argued to 
be a good approximation and yields a simple expression. Convergence is, on the contrary, slower 
in the case of Kg and one would need several multiplets of resonances in order to adequately 
estimate Kg. Fortunately, as far as DT violation is concerned. Kg happens to be multiplied 
by a very small coefficient. Finally, all the EEC’s which are important for DT violation were 
estimated except for the combination K^ + Kq. The difficulty of estimating K^ + Kq comes 
from the complexity of modelling the QCD four-point functions < AAVV > and < VVVV > 
in terms of resonances. One can think of a number of contributions, for instance tensor meson 
exchange, which were never accounted for in DT violation. We have contented ourselves with a 
simple order of magnitude estimate which, after all, could be more reliable than an incomplete 
evaluation. 
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